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Static and dynamic phenomena in overheated free-standing smectic-A films are studied theoret- 
ically. The work is based on a generalization, introduced recently by the authors, of de Gennes' 
theory for a confined presmectic liquid. In this approach, smectic ordering in an overheated film is 
caused by an intrinsic surface contribution to the film free energy and vanishes at some temperature 
depending on the number of layers. Here the theory is further generalized to study the dynamics 
of films with planar inhomogeneities. A static application is to determine the profile of the film 
meniscus and the meniscus contact angle, the results being compared with those of a recent study 
, employing de Gennes' original theory. The dynamical generalization of the theory is based on a 

time-dependent Ginzburg-Landau approach. This is used to compare two modes for layer-thinning 
transitions in overheated free-standing films, namely "uniform thinning" vs. nucleation of disloca- 
tion loops. It is concluded that the nucleation mechanism dominates provided there is a sufficiently 
large pressure difference arising from meniscus curvature. Properties such as the line tension and 
velocity of a moving dislocation line are evaluated self-consistently by the theory. 
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PACS numbers: 6I.30.-v, 61.30. Jf, 64.60.Ht, 64.70.Md, 68.03.Cd, 68.15.+e, 68.60. Dv 



I. INTRODUCTION 



Free-standing films of several smectic-A liquid-crystalline compounds can be heated above the bulk smectic dis- 
ordering temperature without immediately rupturing, and instead are found to undergo successive layer-by-layer 
thinning transitions as the temperature is increased [Jj]-|jJ. The persistence of smectic layering in an overheated thin 
film is usually attributed to enhanced ordering associated with the free surfaces of the film, as is known to occur in 
other contexts Q. There is not yet, however, a clear consensus on the mechanisms by which layer thinning occurs. 
According to one set of theories , thinning takes place when the smectic layer structure in the middle of a film 

vanishes. In an alternative theory [ |12[ |, supported by experimental studies [ p"3| ], layer thinning occurs by spontaneous 
(SJ \ nucleation of dislocation loops prior to the melting of the layer structure in the film interior. This mechanism is not 
■ necessarily unrelated to the former, since a sufficient reduction in the degree of interior smectic ordering is required 
1 ' for it to proceed. 

One of the key experimental observables is the variation of layer-thinning transition temperatures T C (N) with the 
number of film layers N, which is found to be well fit by the power-law relation N cx t~ v , where t — (T C (N) — Tq)/Tq, 
v 0.70 ± 0.10, and To is close to the bulk transition temperature. Alternative mathematical relations |10[ , ^2| , ^4[ and 

S an upper bound Jll| for T C (N) have been derived from the different theories. With appropriate fitting parameters, 
these alternative relations all turn out to agree well with the power-law expression and, thus, are not able to distinguish 
'-^J ■ between the various mechanisms. 

In this paper, we examine further the connections between the different proposed mechanisms of layer-thinning 
transitions. As in several recent works [ ^0|Jr^JT^ |, our analysis is based on de Gennes' |l5| phenomenological Landau 
theory for a "presmectic" film of a fluid exhibiting a second-order bulk A-nematic (N) transition. More precisely, we 
employ a generalization of that theory recently proposed by the present authors UM- One drawback of de Gennes' 
original model stems from attributing surface-enhanced smectic ordering to an cxtcrnal-ficld-like coupling term of 
constant magnitude, which is more appropriate for a film confined between solid walls. This has the consequence that 
a weak degree of smectic ordering in a thin film is predicted to persist up to arbitrarily high temperatures. In order 
to induce layer thinning, recent studies based on this theory have included the effects of a pressure difference AP 
associated with curvature of the meniscus at the film border pCjJi^| . According to the latter studies, AP must be of a 
sufficient magnitude to cause layer thinning, although the layer-thinning transition temperatures are found to depend 
only logarithmically on AP. The modified version of de Gennes' theory proposed by the present authors [ll| utilizes 
a different form of the surface contribution to the free energy (suggested by older theories of wetting [|l6|), which is 
quadratic rather than linear in the surface order parameter and which restores smectic melting at high temperatures 
without requiring the pressure term. 

According to the theory of Ref. pi] , reviewed here in Sec. ||, the free energy per area of an overheated smectic 
film exhibits a discrete sequence of metastable local minima, whose depths decrease with increasing number of layers 
up to some finite maximum N cr depending on temperature. At or slightly below the temperature for which the 
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free-energy well at N cr vanishes, it was assumed in Ref. |Tl[ that the film would spontaneously thin down to a smaller 
thickness. This assumption makes no statement on "how" layer thinning occurs, which is the question addressed in 
this paper. We proceed by generalizing the theory to allow for planar inhomogeneities in the film (Sec, ill) as well 
as the dynamical evolution of the film, using a time-dependent Ginzburg-Landau approach (Sec. |TV| ). In the static 
limit, the theory provides a description of the contact angle between the film and meniscus, which we compare with 
the recent experimental and theoretical studies of Picano et al. JlJ]. Using the dynamical theory, we investigate the 
nucleation of dislocation loops between film domains of different thicknesses and the subsequent growth of the thinner 
region of the film, which is contrasted with the "uniform-thinning" mechanism for the film to achieve a state of lower 
free energy. The results are mapped onto the conventional nucleation picture, in which the activation free-energy 
barrier to nucleation depends on both the difference in well-depths of the homogeneous film regions and the line 
tension E of the dislocation loop. Here we calculate the line tension E, activation energy, and velocity of a growing 
loop self-consistently in the dynamical model, and evaluate their dependence on temperature and number of layers. 
It is found that thinning via nucleation of dislocation loops preempts the uniform-thinning mechanism provided the 
pressure difference AP due to the meniscus is sufficiently large. Conclusions of our findings, and the fact that they 
are restricted to smectic systems with continuous as opposed to first-order bulk transitions, are discussed in Sec. 0. 



II. THEORY OF UNIFORM FILMS 



In this section we review the generalized de Gennes theory in the case of a uniform planar free-standing smectic 
film: further details can be found in Ref. |l]J| . The film is modelled by a thin liquid slab bounded by two parallel 
surfaces located at z = ±L/2, where L is the film thickness. The degree of smectic order in the film is represented by 
the complex order parameter ^(z), where the real part of ^(z) exp(iqoz) describes spatial modulation of the density. 
Here go = 2it/d, with d being the unstressed smectic layer spacing. The order parameter is parameterized as 



^>(z) = V>( z ) ex P[ — i<fi{z)]i 

where ip{z) is the amplitude and 4>(z) = qou(z) is a phase proportional to the layer displacement u(z). 
The Landau free energy per unit area of the film |17| is taken to be 



(1) 
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+ -r s [tP 2 (L/2) + V 2 (~i/2)] - h s fo&(L/2) + 1>(-L/2)] 



(2) 



where C is an elastic constant. This model free energy generalizes that used in Refs. |10 Lj.|lJ,|f5| by including both 
a quartic term (<?/4)?/> 4 with g > in the bulk free-energy density and quadratic contact terms r a ip (±L/2)/2 jl(| in 
the contribution of the surface layers. As in those works, the bulk free-energy density employed in Eq. (||) is strictly 
applicable only to a system exhibiting a second-order bulk smectic disordering transition, which occurs at r = in 
this mean-field theory. We express r — a(T — T*), where T* denotes the bulk mean-field transition temperature [fl9|| . 
Euler-Lagrange equations determining ip{z) and <fi(z) are obtained by functional minimization of Eq. (|^). 

Two alternative versions of this theory for surface-enhanced ordering of smectic free-standing films have been 
investigated. In the first approach jl0],[l2 14 1, following the original work of de Gennes Jig ], such ordering is induced 
solely by the surface field h s , with r s = 0. These studies have also set g = 0. Refining some steps of these earlier 
analyses, which were restricted to the asymptotic limit L £, we find that the equilibrium free energy of this model 
is given by Q 



f(L) = -h s ^{L/2) 



sinh(L/£) 



C [cosh(£/f) - cos(g £)] 



(3) 



where £ = y/C/r is the bulk correlation length. Note that Eq. (Q) diverges on approaching the bulk critical tempera- 
ture, as does the order parameter ip(z) throughout the film, a consequence of neglecting the non-harmonic term gip 4 in 
Eq. (||) . This feature has been mitigated in some versions of the theory |l0|Ji4| , |2l| by the somewhat ad-hoc procedure 
of fixing ip(L/2) instead of h s . It is also worthwhile to note that the minima of Eq. (||) do not occur precisely at the 
values L = Nd, but at slightly compressed values L < Nd of the film thickness, which is due to an interplay between 
finite-size and surface-ordering effects. This feature also holds for the more general model in Eq. (0). 
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In Ref. |ll|, smectic ordering in an overheated film is attributed to a non-zero value r s < 0, with h s = 0. In the 
following we refer to this as the r s -model. In this case, the Euler-Lagrange equations obtained from Eq. (^) always 
admit a trivial solution ip(z) = describing a disordered state of the film, which is the only solution at sufficiently 
high temperature. The numerical solution of those equations, which consider g ^ 0, is described in Ref. |llf| . The 
variation of the equilibrium free energy f(L) typical for this model is depicted in Fig. ^. The free energy exhibits a 
set of wells with centers situated approximately at L = Nd and depths diminishing with N, while f(L) vanishes over 
non-zero ranges of L between the wells. Although not discernable in the figure, the slope of the free energy smoothly 
approaches zero at the limits where f(L) — > 0. In the temperature range T < T s , where T s = T* + r 2 /(Ca), f(Nd) 
for any N is less than some negative threshold depending on T. For temperatures T > T s , wells with non-zero depths 
occur only up to a finite number N cr . This in the case in Fig. [[], where N cr = 11. It was argued in Ref. jll| that the 
temperature at which the free-energy well for L = N cr d disappears, which we will call the "maximum temperature" 
for an A^ cr -layer film, is an upper limit for the layer-thinning transition temperature T c (N cr ). Films of all N < N cr 
can still exist as metastable states and, in principle, thinning could then occur to any one of these states. 

We note the following scaling of the rvmodel free energy, used in Fig. [j] and in subsequent analysis. On expressing 
distances in units of the layer-spacing d, one sees that the free energy in Eq. (||) can be expressed as f(L, r, g, r s , C) = 
(C/d)f(L/d,r,g,f s ), where f = rd 2 /C, g = gd 2 /C, and f s = r s d/C . In principle, the parameter g could also be 
pulled out from / by suitable scaling of the order parameter -0, but we have found it convenient for numerical analysis 
to set this at the (arbitrary) small value g = 0.01 and leave the scaling of / in the form indicated. 



III. NONUNIFORM FILMS: STATICS 



In order to study dislocation loops and associated phenomena in free-standing films, we need to generalize the 
theory of Sec. O to incorporate planar inhomogeneities in a film. This could be carried out by generalizing the free 
energy of Eq. (H) to include horizontal gradients of the order parameter '3/ and nematic director, consistent with de 
Gennes' more general theory of the N-A transition p^ |, as well as of the film thickness L. Full analysis of this theory 
would require minimizing the resulting free energy with respect to ^ and the director. In the remainder of this paper 
we will examine a simplified phenomenological description using an interface-displacement theory |23f| , in which the 
free energy is expanded in gradients of L. To leading order in these gradients, the film free energy F is given by 



F = / d 2 r ± 



f(L) + ±D(V ± L) 2 



(4) 



Here the horizontal direction is represented by r_L, with Cartesian components (x,y), Vj_ = (V^, V y ), L = L(r±) is 
the spatially-varying film thickness, and f(L) is the equilibrium free energy per area of a film of uniform thickness 
L, obtained from Eq. (|^). The coefficient D characterizes the "stiffness" [£3| of the film surfaces, which in principle 
could be derived by analysis of the underlying de Gennes theory. We expect that D ss 7/2, where 7 is the liquid- vapor 
intcrfacial tension. Generally D could vary with L due to changes in the degree of smectic ordering, but such changes 
are expected to be weak and here we will assume that D is constant. 

An additional term may enter Eq. (|^) due to the existence of a positive pressure difference AP = P a i r — Pu qu id 
across the surface of the meniscus surrounding the film. Such a pressure difference produces a shift Afi in the chemical 
potential of the film molecules from their value at coexistence with the vapor phase across a planar interface |^,^5| . 
This leads to the replacement of f(L) in Eq. ^ by f(L), where 

f(L) = f(L) + APL. (5) 

The main effect of the AP term is to shift the depths and to a slight extent the positions of the smectic minima of 
the effective free energy f(L) with respect to those of f(L) : possibly eliminating minima occurring at large L pj. 

The static free energy Eq. (ffl can be used to obtain the shape of the meniscus at the edge of a free-standing film, 
following the analysis of Ref. pj] based on the original de Gennes model. Related considerations have been applied to 
describe the shape of liquid droplets on a solid substrate j2^j27| . The profile of the meniscus is found by minimizing 
Eq. (Q) with respect to L (replacing / by /), subject to the boundary condition that the main part of the film has a 
thickness L = H « Nd. Assuming that L varies only in the cc-direction, the resulting Euler-Lagrange equation is 

The first integral of this equation, with the boundary condition L(~oo) = H, is 
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which has the implicit solution 



'/(J/) 1/2 . (8) 




L(0) 



Here the origin a; — and corresponding value L(0) > H are arbitrary. Note that no boundary conditions need to be 
specified at large positive x, approaching the meniscus border. This reflects the fact that those boundary conditions 
are non-universal and depend on how the film is created However, the description based on the above equations 
is only valid for small Vj_L; deep in the meniscus, higher-order terms in V ±L should enter jL4|p8| . 

Figure || shows the solutions to Eq. (||) for the rvmodel using the same parameters as in Fig. [|. Two cases are 
shown, the solid line corresponding to AP = while the dotted line corresponds to AP = 0.05C/d 2 . (The reduced 
unit for AP is consistent with that for / discussed at the end of Sec. |. Following from this scaling and Eq. (g), the 
horizontal distance x is expressed in units of R sc = d^J Dd/C.) In both cases, the meniscus profile is fairly smooth. 
At low temperatures, when the smectic-A phase is stable in bulk, the meniscus is expected (25) to consist of a set of 
edge dislocations which change the film thickness by steps of height ps d. The absence of such distinct steps in Fig. || 
reflects the fact that here we are in the overheated regime. 

For T < T s , the function f(L) tends to a non-zero value /(oo) with increasing L. The latter represents the 
contribution of surface- induced smectic ordering to twice the interfacial tension 7 of the liquid- vapor interface of a semi- 
infinite liquid. Neglecting small oscillations of /(£) about /(oo) and defining the meniscus slope angle 29 m dL/dx, 
Eq. @j leads to 

92 = ^M°°)-W) + AP(L-H)]. (9) 

On assuming that D = 7/2 and extrapolating the function in Eq. (|^) down to L — H, we obtain for the contact angle 
6 m between the meniscus and film, 

91 = - [/(oo) - f{H)] , (10) 

7 

independent of AP. This result agrees with that derived in Ref. juj; we note that essentially equivalent relations were 
derived some time ago in the case of soap films [^,^9). Note also that Eq. @ predicts growth of 9 with increasing L. 
The meniscus curvature k is found to be 

^ _ dd _ d6 dL _ d6 2 _ AP 
dx dL dx dL 7 

which is just the Laplace law. 

Using the original de Gennes model free energy, Eq. ( ph , one finds that the contact angle diverges as T approaches 
T* from above, whereas the experimental results of Ref. |[L4| indicate regular behavior in this region. This feature was 
alleviated in the calculations of Ref. H], based on the de Gennes model, by fixing tp(L/2) instead of h s . Although 
this eliminates the divergence of the contact angle in the vicinity of the bulk transition point, the resulting model 
predicts a different anomaly, namely that 9 m vanishes as T — > T* for all N. The divergence of the contact angle 
near the bulk second-order transition in the original de Gennes model can be removed by setting g ^ 0. Here we 
present results for 9 m using the rvmodel, although qualitatively similar results are obtained using r s = with h s 7^ 
and g ^ 0. Figure || shows 6^„, in units of C/(d n f), as a function of r oc (T — T*) for various number of layers N, 
using the free energy f(L) depicted in Fig. It is seen that 9 m remains non-zero on approaching the bulk transition 
temperature and increases with decreasing N, in agreement with experiment JWj. We also find that the contact angle 
vanishes above the maximum temperature for a given number of layers. This is in contrast with the model of Ref. 

which yields a small but non-zero contact angle for arbitrarily large temperature and film thickness. However, 
as seen in Fig. ^, the maximum in 9 m as a function of temperature is fairly insensitive to the number of layers, which 
does not accord with the experimental results. 

Comparison of Fig. || with the experimental data G3] for 91 indicates that the scaling unit C/(d r y) should be on 
the order of 10~ 2 . Taking D — 7/2, we then estimate the in-plane distance scale unit to be R sc = d^J Dd/C w 
(10/V2)d fa 2 x 10" 8 m, where we have used the value d = 3 x 10" 9 m M}. 
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IV. DYNAMICS OF NONUNIFORM FILMS 



A. Time-dependent Ginzburg-Landau equation 



To study the dynamics of a thinning free-standing smectic film, we consider the latter to be a two-dimensional 
object with a free energy characterized by its thickness L, given by Eq. (|j) with, generally, / replaced by /. We will 
focus on dynamical processes with large enough characteristic times to neglect inertial effects. Although the details 
of relaxation are undoubtedly quite complicated, here we will proceed by assuming the simplest possible dissipative 
dynamics for L, based on a time-dependent Ginzburg-Landau (TDGL) equation pl| ]. This equation is appropriate 
for describing the dynamics of a non-conserved variable, which L can be regarded in the case of a film open to the 
exchange of molecules with the meniscus. The TDGL equation is 



dL SF 



(12) 



Here 6F/6L is the functional derivative of F, giving the thermodynamic force which drives the system toward equi- 
librium, t is time and 77 is a kinetic coefficient which we will assume to be constant. 

Two opposing mechanisms for thinning of an overheated smectic film are "uniform thinning" (i.e., with V ±L = 0) 
and via nucleation of dislocation loops. If the initial film thickness Lq « Nd is uniform and at a local minimum of 
the shifted free energy /, then it will remain so indefinitely according to Eq. (|l2|). (This picture neglects possible 
disruption due to thermal fluctuations, which we neglect in this work.) Under small displacements of the thickness 
from the initial value Lq, the film will be restored to that initial thickness. Hence uniform thinning can only occur at 
the maximum temperature of an TV-layer film, and then only if AP > ]3^ ]. 

To examine the growth of dislocation loops, for simplicity we first consider a one-dimensional solution of Eq. jl^), 
in the form of an infinite straight-line kink separating domains of thicknesses L\ and Li and moving along the x-axis: 

L = ${x-v 00 t), (13) 



where the function <i> and kink velocity v^, are to be determined. Substituting Eq. (13) into Eq. ( fl2| ) yields the 
ordinary differential equation 



D<f" + V v oo ^-^l = , (14) 



with boundary conditions 



$(-00) =£1, (15a) 
$(00) = Lt, (15b) 

where the prime symbols (/) denote derivatives of $ with respect to its argument. The thicknesses L\ and L2 are 
at local minima of f{L). Usually we will take these to be adjacent minima, with L\ « (N — l)d, L2 ~ Nd, and 
/(L x ) < f(L 2 )_ 

Equation ( |14| ) has a well-known mechanical analogy p3|,p4[ . It can be considered as the dynamical equation 
describing movement of a particle of "mass" D in a medium with "friction coefficient" 77^00 and subject to a "potential 
energy" — f(L). Note that there is no stationary solution (vqq = 0) of this equation unless the depths of the minima 
of f(L) are equal. For any function f{L) characterized by two unequal adjacent minima, as is the case here, there 
should be a unique solution for the velocity Voo > and function $ describing the profile of the kink, which moves 
toward the region of thickness L2 in order to eliminate the domain of higher free energy. Under conditions where the 
higher minimum of f(L) at L2 vanishes and becomes a point of zero curvature, as happens in the present model at 
the maximum temperatures, the velocity Voc and kink shape may become non-unique. In other contexts p3] , |35| , ^6[ , 
this is called a state of marginal stability. Here, we always find (Sec. IV C| ) that uniform thinning occurs under these 
conditions when AP > |3^] . 

The solutions of the one-dimensional equation Eq. (|lj) turn out to be relevant in the more general two-dimensional 
case, as discussed in the next subsection. 
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B. Nucleation of dislocation loops 



If a dislocation loop separating N- and (N — l)-layer re gio ns is nucleated, initially it will be of finite size. According 
to the conventional nucleation picture (see, e.g., Refs. ]l2| , p2| , pgfl for the case of smectic films), the loop then will either 
expand or collapse depending on whether its initial radius is greater or smaller than some "critical" value. Here we 
are interested in determining the critical radius, the associated activation free energy, and the subsequent dynamical 
evolution of the loop. 

We assume that the dislocation loop is a circle. Using in-plane polar coordinates, with origin at the center of the 
loop and radial distance denoted r, Eq. (|l2|) becomes 



dL 



D 



Id ( dL 



r dr \ dr 



df 
dL 



The associated boundary conditions are 



8L{r) 



dr 
L(oo) 



•=o 
L 2 



Nd 



(16) 

(17a) 
(17b) 



The film thickness in the center of the loop at r 
free energy due to formation of the loop is 



will usually be close to the value L\ w (N — l)d. The change in 



AF = ^ / d 2 j 



2/(L) - 2/(L 2 ) + £>(V x i) 



7r / drr 

o 



2f(L)-2f(L 2 )+D 



dL 
dr 



(18) 



In the case of a stationary solution of Eq. (|l6|), corresponding to a "critical" nucleus, AF is the activation free energy. 

The profile of L(r) describing a dislocation loop should have a kink-like structure, with dL/dr w everywhere 
except within a narrow region centered around some value r = R. Numerical solution of the stationary limit of 
Eq. (p^|), with (dL/dt) = and the boundary conditions Eq. ( |l7| ) is difficult, precisely because this is associated 
with a unique "critical" value of the loop radius R. We have found it more expedient, and of more general relevance, 
to solve the full time-dependent partial differential equation Eq. (|l6|), using a standard subroutine (NAG Fortran 
D03PCF). Numerical analysis of the differential equation, starting from initial trial profiles mimicking the expected 
kink structure, shows that solutions in the form of moving kinks are, indeed, dynamically stable. Hence it is worthwhile 
to consider from the outset a solution to Eq. ( |l6| ) in the form of a moving kink: 



L(r,t) =<f>(r-R(t)) 



(19) 



As discussed some time ago in a general context by Chan 
should be a good approximation if the kink radius R is much 
Eq. (Hq) and approximating 1/r by 1/R, we arrive again at Eq 



7|, this form is not an exact solution of Eq. (|1 
arger than its width denoted AR. Using Eq. (11 
Tl), but with Voo defined as |B7|| 




dR 
~dt 



D 1 
~U~R' 



(20) 



This transforms the two-dimensional problem to the one-dimensional case, which requires finding the pair (^foo) 
with constant Vca, satisfying Eq. ( [l4| ) with boundary conditions Eqs. (|l^). The quantity is the asymptotic velocity 
of a loop with sufficiently large radius R and depends (for given D and 77) only on the function f(L) and its chosen 
pair of minima. Our numerical analyses of the original two-dimensional equation Eq. (|l6| ) show that both the kink 
shape $ in the moving coordinate frame and the righ t-han d-side of Eq. ( pp[ ) remain practically constant as the kink 
moves. Thus, while the results to be reported in Sec. IV C are all obtained from the numerical solution of Eq. (|l^), 
the one-dimensional mapping based on Eqs. (0), (fL9|) and (20) usually is an excellent approximation and provides, 



as discussed below, a useful framework for interpreting the results. 
Note that setting the time derivative of R equal to zero in Eq. (E 



Rr = 



D 

VVoc 



j) yields an expression for the critical loop radius, 

(21) 



G 



If Voo is known, then R(t) at arbitrary time can be found by solving Eq. (12G) |37|], which gives the relation 



1 



t - t = / dR' 



R 



R' 



R'-R c 
R - i? + Rc In 



R— R r 



Ra — Rc 



(22) 



where Rq is the radius at an arbitrary initial time to. 

In conventional treatments of nucleation phenomena, properties of the critical nucleus such as R c are expressed in 
terms of purely thermodynamic quantities, and it is worthwhile to show that connection in the present context. From 
Eq. one finds that <& and Woo satisfy 



d 

dX 



1 



-£>($')2 - /($) 



(23) 



where X 



3 t is the argument of <&. Integrating Eq. ( p3| ) over X using the boundary conditions Eq. ([H]) yields 

D 



/(L 2 )-/(L X ) 



where the quantity iJ, which will be identified as the line tension of the dislocation loop, is defined as 

h>\2 



E = D dX(&y = D 



d<f>&. 



(24) 



(25) 



Consistent with the one-dimensional mapping, the lower limit of integration over X has been set to — oo. To a first 
approximation, Eq. ( p4j ) shows that Voa is proportional to the free energy difference fiL^) — f{L\) and, hence, vanishes 
if that difference is zero. Using Eq. (124), Eq. (J2l|) becomes 



Rr = -~ 



E 



f{L 2 ) - /(LO 



(26) 



One can show that Eqs. (|25| ) and ( pq ) for E and R c are consistent with standard arguments of nucleation theory. The 
free energy change Eq. (|18|), with the solution for L in the form of Eq. ([l^), can be approximated as 



AF = -irR 2 [f(L 2 ) - f{L x )] + 2ttRE 



(27) 



This follows on approximating the integrand in Eq. (|l8| ) by the constant f(L\) — f(L 2 ) inside a circle of radius 
R — A_R/2, neglecting the integrand for r > R + AR/2, and using Eq. (|2^) within the kink region of width AR. 
Equation (|2^), with E given by Eq. (p5[) , is then obtained on assuming AR <C fi- The first term in Eq. (|27]) is the 
decrease in film free energy due to the difference /(-f^) — |58| while the second term is the free energy which 

has to be overcome due to the line tension of the loop. Maximizing AF with respect to R yields the critical radius 
R c given by Eq. (p(|). The corresponding activation free energy F act = AF(R C ) is 



Fact 



ttE 2 



f(L 2 ) - /(LO 



itR c E. 



(28) 



The expression for E in Ec 



250 agrees with a familar mean-field relation for the tension of a stationary interface 



in terms of its profile shape [19|. Here that relation also applies to a moving kink of sufficiently large radius, under 
the assumption (supported by our numerical studies) that the profile shape is preserved during its motion. 

One final point to note concerns the physical interpretation of the equation of motion Eq. (^) for R(t). That 
equation is equivalent to the balance of thermodynamic and dissipative forces per unit length of the dislocation line: 



I / dAF \ 
2^R V dR J 



E (dR\ 



(29) 



This yields Eq. (|20|) on using Eqs. (g4J) and ([27|), and agrees with the model of dislocation-loop dynamics described 
by Geminard et al. Pq| , on identifying the mobility /i used by these authors with the quantity Dd/(r]E). The study 
of dislocation-loop dynamics in Ref. |28[ |, performed below the bulk A-N transition temperature, showed that the 
resulting integrated form Eq. ( p2| ) very well fits experimental data. Performing our own fit of Eq. ( p2| ) to the data 
reported in Ref. |^8| yields the values Vqo = 2.59/im/s and R c — 42.6/mi, to be contrasted with results discussed below 
for overheated films. 
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C. Numerical results 



Figure ^ illustrates the dependence of the free energy AF on radius R due to a growing dislocation loop in a 6-layer 
film, for the case of /(£) shown in Fig. [j] and AP = 0. The in-plane distance is expressed in units of R sc = dy/ Dd/C, 



as in Fig. |L while the free energy is plotted in units of Dd 2 , which follows by scaling of Eq. (18|). The crosses are 
obtained by monitoring the numerical solution of Eq. ( |l6| ) for a dynamically stable kink as it evolves with time, where 
the loop radius R is defined such that L(R) = [L(0) + L(oo)]/2. It is seen that AF(R) is very well fit by Eq. j27|): 
the fitted value of /(.L2) — agrees with that obtained directly from the static theory with a precision of 0.5 

per cent. We also verified that the right-hand-side of Eq. (20) remains constant within the same precision, giving 
a value Voo = 0.18D/(r)R sc ). This numerical calculation supports the ansatz E q. fp jjh for a two-dimensional kink 
with constant Voo related to R(t) by Eq. (^0|). Analogous fitting to Eqs. ( f27|) and (|20|) of the numerically determined 
AF(R) and kink radius has been carried out for all our calculations, showing comparable precision except very close 
to the maximum temperature of an iV-layer film, where uniform thinning occurs. As a further check, we find very 
close agreement between the values of the line tension E obtained from the fitting to Eq. ( p7| ) and by direct evaluation 
of Eq. (|2^) from the numerically determined kink profiles. 

Only recently, in Ref. |l3j , have the dynamics of dislocation loops in overheated films been studied experimentally, 
albeit for a system exhibiting a first- order bulk A-isotropic (I) transition. (This difference will be discussed in Sec. ^.) 
The magnitude of the loop velocity was found to be 10 3 -10 4 times larger than that reported for the smectic-A phase 
in Ref. p8| , which we will comment on further below. We note that the data in Ref. jl3| show a purely linear 
dependence of R on time, suggesting that R/R c » 1 in the experimentally accessible range and that the measured 
velocity corresponds to of the present theory. The data reported in Ref. |l3| also show that the dislocation-loop 
velocity slightly increases with increasing number of layers. To check the dependence of Voo on thickness in our theory 
we took r] to be independent of N. Figure || shows the results of calculations using the f(L) function employed in 
Fig. 0. It is seen that in some range of thickness the velocity slightly increases with N provided that AP > 0. We 
do not rule out that rj depends on thickness and diminishes with N, as suggested in Ref. but a detailed analysis 
of rj is beyond the scope of this paper. Figure ^| shows, in the case of a 5-layer film, that is predicted to increase 
with temperature, in agreement with experiment |l3|. The rate of increase is enhanced in the presence of a nonzero 
AP. The upper temperature limits of the curves in Fig. [}] are slightly less than the maximum temperatures for the 
given values of N and AP. Beyond these limits, dislocation-loop growth is superceded by uniform thinning. 

Using Eq. (^l|) and the velocity data in Fig. ||, the critical radius R c can be determined. Referring only to the points 
for N = 11 in the figure, we find that R c = U.9R SC and 1.43P SC for AP = and AP = 0.05C/d 2 , respectively. These 
values bracket the range of R c values obtained for other values of N. Using the estimate for R sc described at the end 



of Sec. Ill, we thus find R c to be in the range 10 1 to 10 2 /im, several orders of magnitude smaller than the value 



deduced above in the bulk A phase. In view of Eq. (|21j), these results are consistent with the reported differences in 
loop velocity i»oo below and above the bulk transition temperature. 

The key quantity which determines whether spontaneous nucleation of a dislocation loop in an overheated film 
actually occurs is the activation free energy F act . The dependence of F act (in units of Dd 2 ), for a single-layer 
dislocation loop, on the number of layers N in the initial uniform film is shown in Fig. |^. Again, the calculations 
were done using the f(L) function depicted in Fig. |l|. Note that the temperature in this case is slightly below the 
maximum temperature for an 11-layer film. To judge whether nucleation occurs, we use the argument described in 
Ref. based on the frequency per unit area of forming a dislocation l oop of radius R > R c . This frequency is given 
by / = fo ex p( — F act /kBT), where /o is estimated to be 10 26 cm~ 2 s _1 (22). This yields |l2| that the condition for a 
dislocation loop to nucleate in 0.1 s in a 1 cm 2 film is F act /kBT m 60. To compare this number with Fig. [7] requires 
a value for the dimensionless parameter Dd 2 /(UbT). Taking D = 7/2 w 1.5 x 10~ 2 N/m, fcsT = 4.5 x 10 _21 J and 
d = 3x 10~ 9 m, we estimate that Dd 2 jkuT w 30. These numbers yield the criterion that F act / (Dd 2 ) must become 
about 2 or smaller for spontaneous nucleation to occur. Figure ^ shows that, for the chosen set of parameters, this 
can be achieved provided that AP > is sufficiently large. 

Figure || shows the variation of the line tension E with initial number of layers N for the same cases shown in 
Figs. [| and [?]. Except for the last two points in the figure (N = 10, 11), E is seen to be essentially independent of 
AP. It follows from Eq. ( p8| ) that it is the decreasing magnitude of E on approaching N cr for a given T which is 
mainly responsible for the decrease in F act in Fig. |?] near N cr and, similarly, for the increase of Voa when AP > 
in Fig. ||. The difference between the two cases, AP = and AP > 0, in Figs. |5| and ^ is due to the fact that 
/(i2) — f(Li) « 7(^2) — f(Li) + APd is dominated by the term APd for the value of the latter used. 
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V. CONCLUSIONS 



The present theory, a modification of de Gennes' fig] theory of presmectic films, is based on the generally accepted 
view [T^ ] that the occurrence of overheated free-standing films is due to surface-enhanced smectic ordering H. In 
the case of a uniform planar film, our theory fll|| predicts that there is a maximum temperature for which smectic 
ordering in an Allayer film can occur. We associate this with an upper bound for the true layer-thinning transition 
temperature T C (N). Employing a dynamical generalization of the theory based on a time-dependent Ginzburg-Landau 
(TDGL) equation, we have shown that thinning via nucleation of dislocation loops, the mechanism indicated by recent 
experiments [[L^-[l4f , is possible provided the pressure difference AP resulting from curvature of the surrounding 
meniscus is sufficiently large. Otherwise, the film would undergo either "uniform thinning" or, possibly, rupturing by 
a process analogous to spinodal dewetting |52| . The requirement for a non-zero AP to promote nucleation of dislocation 
loops is consistent with other recent studies |I(|[l4j], although we emphasize here that the condition AP 7^ is not 
essential for some type of thinning process to occur. 

In the present work, all nucleation properties (i.e., R Cl F act , Woo, E) arc interrelated within the framework of solutions 
to the TDGL equation. In particular, the line tension E of the dislocation loop is expressed in terms of the kink 
profile, see Eq. (|25|). From this equation, E depends implicitly on the elastic behavior of the system through the 
dependence of the kink shape on the uniform film free energy f(L). This self-consistent method for incorporating the 
line tension differs from those used in Refs. flC||I1^ , |l4| , pH . 

It is important to emphasize that the present work is restricted to smectic-A liquids undergoing continuous A-N 
transitions in bulk. For this reason we do not attempt to make quantitative comparisons with the experimental results 
of Ref. [ |l3| for dislocation loop dynamics, since the latter pertain to a system with a first-order A-I transition. Such 
a system can be treated using the present Landau-de Gennes theory (albeit with considerably greater complexity 
ppf ) by appropriately modifying the bulk free-energy density in Eq. (g). One significant difference which is expected 
concerns contributions to the "AP" term in the resulting uniform-film free energy. In a system with a first-order 
bulk transition, such a term arises even in the absence of meniscus effects, due to the grand-canonical [ fL7j free-energy 
difference between a metastable bulk smectic-A phase and the isotropic phase ||[l8|. This effect was recognized in 
the nucleation theory of Ref. |Q, although the latter work otherwise employed de Gennes' theory for a second-order 
bulk transition. Clearly, to compare the present dynamical predictions with experiment, it would be of interest to 
perform measurements on the dynamics of layer thinning in systems such as that studied in Ref. H] , which exhibit 
second-order bulk A-N transitions. 

Our picture is that dislocation-mediated thinning of an overheated free-standing smectic film may preempt the 
uniform-thinning mechanism and thus, for given N and fixed values of the model parameters such as r s , occur at 
a lower temperature than predicted |ll| by considering a purely uniform film. This is conceivable because all free- 
standing film states are metastable |]I8|. Here we have presented several qualitative results in support of this picture. 
An additional step would be to evaluate the shift in layer-thinning transition temperatures T C (N) from those predicted 
jTIJ for a uniform film, and attempt a refitting with experimental data. This is a somewhat tedious task and is left 
for future study |2C|| . Further studies will also be directed to evaluating the kinetic coefficient r\ and, as already 
mentioned, extending the present theory to smectic films with first-order bulk transitions. 
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FIG. 2. Thickness of the meniscus vs. scaled in-plane distance x/R sc , using the f(L) function of Fig. |l} The solid line 
corresponds to AP — while the dotted line corresponds to AP — 0.05C/d 2 . The origin x = is such that, for both curves, 
L(x = 0) = 6.00001 
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FIG. 3. Squared menisus contact angle 8^ in units of C/'yd as a function of dimensionless temperature variable 
rd 2 /C oc (T — T*) for various number of layers N, calculated in the r s -model with the same parameters as in Fig. hi. 
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FIG. 4. Free energy AF (reduced by Dd 2 ) of a dislocation loop in a 6-layer film vs. radius R/R sc Crosses: numerical 
solution of the dynamical equation Eq. jl^ ) for the same f(L) and r s -model parameters as in Fig. 0, with AP = 0. Solid line: 
fit of AF with the parabolic function in Eq. @. The best fit yields E = 0.54V DCd and f(L 2 ) - f{Li) = 0.098C/d. 
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FIG. 5. Velocity i;^ (reduced by D/(-qR sc )) for iV-layer films vs. the number of layers, calculated using the f(L) depicted 
in Fig. 0. Circles: AP = 0. Triangles: AP = 0.05C/d 2 (as in Fig. [|. 
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FIG. 6. Reduced velocity for a 5-layer film vs. rd 2 /C oc (T — T*). Parameters of the r s -model are the same as in Fig 
Lower curve (solid): AP = 0. Upper curve (dashed): AP = 0.05C/d 2 . 
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FIG. 7. Activation free energy F act (reduced by Dd 2 ) for TV-layer films vs. the number of layers calculated using the f{L) 
function shown in Fig. @. Circles: AP = 0. Triangles: AP = 0.05C/d 2 . 
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FIG. 8. Line tension E (reduced by V DCd) for iV-layer films vs. the number of layers, calculated using f(L) depicted in 
Fig. §. Circles: AP = 0. Triangles: AP = 0.05C/d 2 . 
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